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•Abstract 

During the Apollo program two successful heat flow measurements 
were made in situ on the lunar surface. At the Apollo 15 site a value of 

x 10 watts /cm was measured and at the Apollo 17 site a value of 

- 6 2 

2.2 x 10 watts/cm was determined. Both measurements have uncer- 
tainty limits of ± 20% and have been corrected for perturbing topographic 
effects. The apparent difference between the observations may correlate 
with observed variations in the surface abundance of thorium. Comparison 
with earlier .determinations of heat flow, using the microwave emission 
spectrum from the moon, gives support to the high gradients and heat 
flows observed in situ. 
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IN -SITU MEASUREMENTS OF LUNAR HEAT FLOW 


Introduction 

There has long been an interest in the rate at which heat is 
escaping from the moon. The moon is a planetary-sized body that re- 
presents a significant sample of the solar system in the region of the 
ter re stria L planets . On the other hand, the moon is a small enough body 
so that there is reason to believe that during its 4. 6 billion-year history 
it has lost a significant portion of its initial heat, and, as a consequence, 

the present heat flux mainly results from heat generated by radio iso- 

I: • 

topes in the interior to a depth of about 300 km. 

Petrological and geochemical evidence gained from surface samples 
indicates that the moon was radially differentiated early in its history. 

During this differentiation the long-lived, heat-generating isotopes of 

238 23 5 40 232 t 

U,. U,- K and Th were purged from the interior and concen- 

trated in the outer layer of the moon. As a result, the surface heat flow 
from the moon should very nearly reflect the total abundance of these* 
isotopes in the moon, and thereby provide a valuable chemical constraint 
on the moon’s bulk composition. 

Prior to the Apollo missions, lunar heat-flow determinations were 
based on earth -based observations of thermal emissions from the moon 
in the microwave band. Because of the partial transparency of lunar 


material to electromagnetic waves longer than 1 mm, the emission 



spectrum at wavelengths greater than 1 mm depends on temperatures in 
the subsurface* If the electrical properties of lunar soil are known, the 
subsurface temperature profile can be determined from the emission 
spectrum. 

The most comprehensive effort to detect heat flow from the interior 

by this technique has been made by Troitsky and colleagues (1, 2) at the 

Badiophysical Research Institute, Gorky, U. S.S. B. Their well-known 

curve, shown in Figure 1, indicates an increase in brightness temperature 

with wavelength of about 0. 6 C/cm, Using electrical and thermal 

properties deduced from microwave observations in the 1mm to 3cm 

range, Tikhonova and Troitsky (2) interpreted this spectral gradient in 

of 

terms of a heat flow/3 x 10 to 4 x 10 ^ W/cm^. Such a heat flow is 
approximately 1/2 the mean of observed heat flow values On the earth. 

In Situ Measurements During the Apollo Program; 

The manned lunar landings of the Apollo program provided an 
opportunity to make direct measurements in the lunar surface layer 
relevant to the heat flow through the surface. Successful measurements 
were made at two of the landing sites: Hadley Bille, near the edge of 
the Imbrium basin - visited on Apollo 15 - and Taurus Littrow, a narrow 

embayment on the southeastern margin of Serenitatis - visited on Apollo 
17. See Figure 2. 

At each location the astronauts buried two probes in the lunar soil 
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to measure the temperature and thermal conductivity of the soil. At 
the Apollo 15 site the probes were buried to depths of 1. 0 and 1.4 m. 
and at Apollo 17 both probes wpre buried to a depth of 2. 3 m. Each 
probe contains eight platinum resistance thermometers and four thermo- 
couples which detect temperature at 11 different levels in the subsurface. 
Four thermometers on each probe are surrounded by heaters which can 
be turned on by command from earth. These heaters are used to make 
in situ determinations of thermal conductivity. The range and accuracy 
of measurements made by the heat flow experiment. are shown in Table 1 
The platinum resistance thermometers were carefully tested to demons- 
trate that they would retain their calibrations after experiencing the 

mechanical and thermal shocks of the lunar mission. Temperature 
data from all the thermometers are relayed to earth every 7. 2 minutes. 
At the Apollo 15 site we presently have more than two: and a half years 
of data and more than a year of data at the Apollo 17 site. 

A Summary of Results 

The experiments installed on- the moon provide extensive infor- 
mation on the temperature and thermal properties of the lunar surface 
layer, to a depth of; 3 meters, including surface temperature variations , 
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near-surface thermal properties, subsurface temperature variations 
and thermal conductivity. All of this information is essential to under- 
stand the total heat budget near the lunar surface and the contribution 
of the flux from the interior. 

Su rface Temperature Variations ; Measurements by thermocouples in 
cables above the lunar surface provide information on the surface tem- 
perature variation. The cable is in radiative equilibrium with the lunar 
surface, except during times when the temperature is changing rapidly 
as during an eclipse or at terminator crossing. The lunar surface tem- 
perature can be readily computed from laws governing thermal radiation. 

In Figure 3 we show the surface temperature variation at the Apollo 17 
site during a complete lunation. During lunar day the temperature 
deductions have large errors because of uncertainties in the amount of 
solar radiation reflected from the lunar surface, but at night the errors 
are small. Similarly, surface temperatures can be deduced quite accurately 
from thermocouple data during the umbral phase of an eclipse. 

, In a manner similar to the classical methods of Wesselink (3) and 
Jaeger (4), the cool-down of the surface after sundown and during an 
eclipse can be used to deduce the thermal properties of the regolith to a 
depth of about 15 cm. For analysis of the in situ data we use a thermal 
model that includes many layers with thermal properties that vary with 

1 

depth and temperature. To explain the observed temperature variations 
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at. lunar night and during an eclipse, the conductivity and density must 
vary with depth. The variations of density and conductivity shown in 
Figure 4 will explain the surface temperature variations during the 
lunar night almost exactly but these deduced profiles are not necessarily 
unique. Two features of the profiles shown are essential to explain the 
data 

1. . The upper 1 to 2 cm must have an extremely low thermal conductivity 
and this conductivity must be temperature dependent. The conductivity 

at the mean surface temperature (216 fC) is approximately 1.5x10 W/cm-K, 
which is in good agreement with measurements on returned lunar fines. 

2. At a depth of about 2 cm the conductivity must increase greatly to 
values 5 to 7 times greater than the surface value. 

The Near-Surface Mean Temperature Gradient: One of the most interesting 

features of the subsurface temperature measurements is the very large 

difference in mean temperature (i. e. the temperature averaged over one 

centimeters 

lunation) between the surface and depths, of a few / ; . At the Apollo 15 
site the mean temperature 35 cm below the surface is 45°K. higher than 
at the surface and the difference at the Apollo 17 site is 40°K. This large 
increase in mean temperature is due primarily to the temperature 
dependence of thermal conductivity in the top 1 to 2 cm, which results from 
the non-linear behavior of radiative heat transfer. These large differences 

' V? t , 

require that the ratio of the radiative component to conductive component 
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at 35>0°K must be about 2 to 3. Figure 5, from Keihm and Langseth (5) 
shows the variation of mean temperature and the amplitude and phase of 
variations of lunation period with depth in the regolith based on the models 
shown in Figure 4. 

Subsurface Temperature Profiles 

The probes are inserted inside the hollow fiberglass tubing which 
is drilled into the lunar soil. Figure 6 shows Charles Duke, an Apollo 16 
astronaut drilling one of the heat flow holes. Figure 7 shows the tem- 
perature history of one of the probes after insertion into the tube. The 
probes require approximately one month to reach within a few thousandths 
of a degree of thermal equilibrium with the surrounding regolith. The 
thermometers buried below 80 cm depths do not see any perceptible 
variation due to the monthly temperature cycle, and temperature gradients 
should reflect heat flowing from the lunar crust. The temperature profiles 
at four of the probes are shown in Figure 8. Very small corrections have 
been added to these data to account for thermal shunting effects of the fiber- 
glass tubes and probes so that: these profiles should represent true undis- 
turbed regolith temperatures. 

Thermal Conductivity: Thermal conductivity of the regolith can be deduced 

from three different effects. First, and most important, were measurements 
made by in situ experiments. The effects of heaters turned on for a period 
of 36 hours at low power were measured. From the rate of rise of tem- 
perature after 20 hours it is possible to determine the conductivity. 
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Second, the initial cool-down of the probes from high temperatures 
permit determination of conductivity based on the initial energy input 
into the hole. Cool-down estimates can be made at all the gradient 
sensors, i. e. at eight different depths in each hole. Third, tempera- 
ture variations with a monthly period penetrate to approximately 80 cm 
and the annual variation of surface temperature is felt at all depths. 

The attenuation of these variations with depth depends in part on the thermal 
conductivity of the surrounding material. However, because of radiative 
transfer along the fiberglass tubing, the attenuation data are difficult to 
interpret. Our analysis requires a thermal conductivity between 1 
and 2 x 10 ^watts/cm-°K to reproduce the attenuation observed. The con- 
ductivities measured by the first two techniques are shown versus depth 
in Figure 9. We note that the thermal conductivity of the lunar soil lies 
between 1. 4 and 3. 0 x 10 watt/cm- K;. this is approximately a factor 
of 10 higher than the conductivity at the surface. The increase of con- 
ductivity at about 2 cm depth appears to be mainly due to a large increase 
in the soil compaction and grain boundary contacts with depth. It is 
likely that at this 2 cm depth the disruptive effects of micrometeorite 

bombardment give way to compactive effects. Conductivity values for 

-4 o 

regolith fines as high as 2 x 10 watts /cm- K have not been duplicated 
in the laboratory and further tests of highly compacted lunar soil should 
be made. 


t 


Heat Flow In Situ: When the conductivity and gradient observations are 


combined, the heat-flow values shown in Table 2 result. The best 

-6 2 

value of heat flow at the Hadley Rille site is 3. 1 x 10 watts /cm and 
-6 2 

2. 8 x 10 watts /cm at Taurus Littrow. 

Temperature measurements at probe #2 at the Apollo 17 site 

require special attention, see Figure 8. The profile indicates a very 

large decrease in gradient with depth at 130 cm which, because the 

conductivity is relatively uniform, must reflect a change in heat flow 

with depth. Also, the heat flow in the lower part of the hole, about 

1,9 watts /cm^, is considerably lower than that at the other two locations. 

These results suggest that heat flow is locally disturbed, perhaps by 

a large rock buried very near where the probe is emplaced. The heat 

flow, using the temperatures at 67 and 234 cm, gives a heat flow of 2. 5 
- 6 2 

x 10 watts /cm , which is in reasonable agreement with the value at 

. • 

probe #1. 

Possible Disturbance to the Heat Flow: 

How representative are the measurements of the average heat 
loss from the moon? The answer to this question depends on whether 
there are significant regional and local disturbances. Certain disturbing 
effects, such as that of local topography, can be estimated and corrected 

The amounts of radioisotopes in the crust may be anomalous in the 


for. 
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region where the heat-flow observations are made. We have orbital 
data on the distribution of thorium and uranium on the surface which 
can be applied to this problem. 

Other effects, such as thermal conductivity contrasts in the sub- 
surface which can refract the heat-flow lines, are not directly obser- 
vable. but geological observations can be used as a guide for assuming 
subsurface conductivity geometries, and from these assumptions an 
appreciation of whether the measurements are anomalous may be 

obtained. 

Topography : Significant disturbances of heat flow will occur m the 

vicinity of craters which have a diameter-to-depth ratio of 6 or less. 

The dominant effect of such craters is to increase the heat flow just 
outside the crater rim due to the slightly higher mean temperature in 
the crater floor. In the upper part of Figure 10 we show the heat-flow 
anomaly over a crater with a diameter-to-depth ratio of 6 as a function 
of radius. It can be seen that the anomaly decreases very rapidly with 
distance from the rim. By and large the astronauts were successful in 
setting up the experiment far from craters larger than 1 m in diameter. 

At the Apollo 17 site topographic maps show three craters which we 
estimate have a combined effect that increases the heat flow by 0. 3 W/cm' 
That is, a correction of about -10% should be applied to the Taurus 
Littrow values for the effect of craters. At the Hadley Rille site there 



were no craters in the vicinity of the probes which would have a signi- 


ficant effect on the heat flow. Consequently, no correction has been 
applied. 

At Hadley Rille both the rille and the Apennine Front will affect 

the heat flow, but in opposite ways. Both effects are on the order of 5% 

and thus appear to be self cancelling. Thus, it appears that the best 

value for the heat flow at Hadley Rille is the uncorrected value which is 
-6 2 

3, 1 x 10 watts /cm , with an estimated uncertainty of :fc 20%. 

The massifs that bound the Taurus Littrow valley on the north and 

south have a significant effect on the heat flow. We have estimated the 

correction to be applied using a method developed by Lachenbruch (6). 

The valley is modeled as shown at the bottom of Figure 10. Based on 

this model we estimate that a correction of -15% to -20% should be 

applied to the Apollo 17 measurement. Applying all corrections, the 

the region ^ . » 

best value for heat flow in/Taurus Littrow/is 2. 2x 10 watts/cm with 

an estimated uncertainty of ± 20%. 

Surface Radioactivity : When topographic effects are taken into account 

the heat flow at Taurus Littrow is 25% to 30% lower than at Hadley Rille. 
The results of the orbiting gamma ray experiment reported by Metzger 
(7) gave evidence of substantial variations in radioactive elements on 
the surface. One of the regions with the highest concentrations in radio- 
activity is the Hadley Rille area of Mare Imbrium. There the counts 



per second are about twice those at Taurus Littrow and 3 to 4 times 
those observed over much of the lunar farside.. The difference in heat 
flow between these two sites may therefore reflect a real variation in 
radioactive heat production in the lunar crust, A similar correlation 
between surface heat flow and radioactive heat production of surface 
rocks is observed on the earth. A most significant implication is that 
two in situ measurements may overestimate global heat flow, especially 
if the results of the gamma ray experiment over the farside and nearside 
highlands are representative since the observed concentrations of radio- 
active isotopes are much lower there. 

Possible Subsurface Effects: Both sites where the heat flow experiment 

were installed are at the margins of large raascon basins. These basins 
have been flooded by basaltic lavas early in the moon's history. and as a 
result it is possible that a conductivity contrast exists between the mare 
basalt and the underlying basin floor and adjacent highlands. It is not 
likely that the bulk conductivity of the basalt will be as high as that 
measured on returned solid samples. The active seismic experiments by 
Kovach and Watkins (8) indicate that these flows are highly fractured. 
Extensive fracturing of rocks in vacuum will decrease their conductivity 
appreciably. It may be possible to entertain a conductivity contrast of 
mare fill material to underlying basin floor material of 10. 

Our observations which are located at the margin of basins may see 


an edge effect as a consequence of this contrast. If they lie within the 
basin as the Apollo 15 site appears to, an edge effect heat flow anomaly 
would add to the observed heat flow. If on the other hand the observation 
is outside the basin rim as might be the case at Taurus Littrow, the 
observed heat flow would be disturbed toward lower values. Uncertainties 
do not at this time permit an accurate assessment of disturbance due to 
buried contacts between rocks of different thermal conductivity. The 
best estimate of the size of such an effect is by comparison with other 
estimates of heat flow from the moon, for example, those made from 
earth-based microwave measurements. 

Comparison With Earth-Based Microwave Measurements 

To compare our measurements with microwave measurements we 
return to the set of measurements of lunar brightness temperature between 
3 and 70 cm made by Troitsky and colleagues. Waves from 5 to 20 cm 
are emitted from depths comparable to those measured by the heat-flow 
experiment. We will compare the spectral gradient in this band of wave- 
lengths with that expected f^orq a lunar surface layer with the temperatures 
and thermal properties we measured in situ. The greatest uncertainty 
in such a comparison is the value of the power absorption length, 
of electromagnetic waves. In Figure 11 we show the microwave data 
compared with theoretical results using different absorption lengths. The 
model has the parameters given in the Table at the top of the figure. 



The lowermost curve has a power absorption length which is a function 
' . ^ '1.48 

of \ i. e. £ e " 5. 8 a which is an empirical fit to the experimental 
results. This relation produces a good fit to radiotelcscopic obser- 
vations of the attenuation of the variations of microwaves over. a monthly 
period in the range from 1mm to 3. 2cm. The result is an increase in 

brightness temperature that fits within the error bars but is a poor fit 
spectral gradient 

to the observed/ . The temperature gradient in the lunar surface layer 

would have to be larger to better fit the microwave data using this relation. 

The middle curve corresponds to an of 50 X , which is near the mean 

of values reported by Gold et at. (9) based on measurements of lunar 

.spectral 

samples ata wavelength of 68 cm. The best fit to the/gradient would be 
obtained for an of 80)s . The principal result of this comparison at 
this point in our knowledge is that gradients of 1. 3°C /m or higher as 
observed in situ are supported by microwave observations. The heat 

< ' ’ ?• ' -A 7 

flows deduced by Tikhonova and Troitsky from 2. 9 x 10.” watts /cm to 
-6 2 

4x 10 watts/cm are in close agreement with our results. 

Future Work; -• 

Several lines of future work can be proposed. 

First, more laboratory measurements of the electromagnetic 
of lunar fines 

power absorption lengtlyin the wave band from 5 to 20 cm 
would greatly improve the comparisons made above. 

Second, further earth-based measurements of microwave emissions 



in the wavelength range from 5 to 30 cm would be extremely valuable, 
Measurements with sufficient resolution to detect variations in emission 
spectra over the lunar disk would be especially significant. 

Third, other in situ measurements in highland regions, using an 
automated lander, would be extremely important. 

Last, microwave observations from lunar orbit could possibly 
map variation in heat flow over the whole moon. 
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TABLE 1 


HEAT FLOW EXPERIMENT TEMPERATURE MEASUREMENTS 


MEASUREMENT 

RANGE 

ACCURACY 

Platinum Resistance Thermometers: 
Absolute Temperature 

190-270°K 

±0: 05°K 

Temperature Difference 

±2°K 

±0. 001°K 

Cable Thermocouples; 

Absolute Temperature- * . 

70~400°K 

±0. 5° K 













SUMMARY OF LUNAR HEAT FLOW RESULTS 



Gradient 

°K/m 

Conductivity 

WxlO _4 /cm°K 

Heat Flow 

WxUrfycm 2 

Apollo 15 
Probe 1 

1.75 

1.78 

3. 1 

Apollo 17 
Probe 1 

1.3 6 

2. 06 

2. 8 

Probe 2 

1. 30 

2. 00 

Z. 5 

1 T -jAfirf. 


The estimated error of heat flow determinations is ±20% 




Figure 1 


Figure 2 


Figure 3 


Figure 4 
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Figure 5 
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FIGURE CAPTIONS 

Data from microwave observations made in the 
Soviet Union between 1961-1964 using the "arti- 
ficial moon" technique (1,2). The plot shows 
microwave brightness temperature versus 
wavelength of the radiation. 

A lunar map showing the locations of the two 
successful heat flow measurements. 

1 ' 

The calculated surface temperature throughout a 
lunation based on the temperature of a thermo- 
couple inside a cable (see inset) which is exposed 
above the surface. 

4 f ’ 

Profiles of thermal conductivity and density with 

depth in the lunar regolith which will explain the 

observed surface temperature variations shown 

in Figure 3. Conductivity values shown are those 
appropriate to the mean temperature at each depth. 

The plot to the right shows the expected peak to 

peak monthly variation of temperature as a function 


t 
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Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


of depth. The left hand plot shows mean tem- 
perature and phase lag of the monthly variation 
versus depth. These results are based on the 
conductivity and density models shown in Figure 4. 

An astronaut is shown drilling sections of fiber- 
glass tubing into the moon for the heat flow probes. 

The 10 month temperature history of eight thermo- 
meters on one of the probes at the Apollo 17 site. 
The inset shows the temperature depth profile 
after 75 days. 

Temperature depth curves for the four lunar heat 

flow measurements. The hatched areas above 

70 cm are the envelopes of monthly variations. 

* ♦ 

A summary of thermal conductivity determinations 
at the heat flow site. At the top of the plot the 
lines are the same as in Figure 4. Deeper points 
are results of in situ measurements and analysis 
of probe cool-down data. 

In the upper part of the figure the heat flow anomaly 
expected over a crater with an aspect ratio of 1. 6 


is shown. The lower part shows schematically 
a cross section of Taurus Littrow and the estimated 
effect on the heat flow, Aq. 

Figure 11 The data shown in Figure 1 compared with theo- 

retical emission spectra based on the thermal 
properties models at the Apollo 15 site. Three 
different relations between absorption length 
. and wavelength were used and are noted on the 


curves. 
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